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Abstract 

We investigate the fermionic decays of top squarks ii^2 and bottom squarks 61^2 
in the Minimal Supersymmetric Standard Model with complex parameters Mi, fi, At 
and Af) . In the analysis we particularly take into account the cosmological bounds 
imposed by WMAP data. We plot the CP phase dependences of stop and sbottom 
decay widths. 

PACS numbers: 14.80.Ly, 12.60. Jv 

1 Introduction 

SUSY provides an appealing realization of the Higgs mechanism for mass generation. When 
SUSY is broken softly, superpartners acquire masses not exceeding 1 TeV. Since these new 
particles are in the exploration range of the LHC, we have to analyse non-standard Higgs 
sector extensively in all possible ways [1]. 

The minimal supersymmetric extension of the Standard Model (MSSM) requires a non- 
minimal Higgs sector [2] which introduces additional sources of CP-violation beyond the 
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economical Kobayashi-Maskawa phase of the SM. The plethora of CP- phases also influences 
the decays and mixings of B mesons (as well as D and K mesons). The present experiments 
at BABAR, Tevatron and KEK and the one to start at the LHC will be able to measure 
various decay channels to determine if there are supersymmetric sources of CP violation. 
In particular, CP-asymmetry and decay rate of S — > form a good testing ground for 
low-energy supersymmetry with CP violation [3] . These additional sources of CP- violation 
are welcome to explain the cosmological baryon asymmetry of the universe. In addition 
to this, the lightest superpartner, i.e. the lightest neutrahno Xi could be an excellent 
candidate for cold dark matter in the universe. 

With the precision experiments by Wilkinson Microwave Anisotropy Probe (WMAP), 
the relic density of cold dark matter can be constrained to 0.0945 < flcDMh^ < 0.1287 at 
2(7 level [4]. Recently, in the light of this cosmological constraint, an extensive analysis of 
the neutralino relic density in the presence of CP phases has been given by Belanger et 
al. [5]. 

Analyses of the decays of third generation scalar quarks with complex SUSY parameters 
have been performed by Bartl et al. [6]. In the present note we repeat this analysis taking 
into account the cosmological bound imposed by WMAP. Namely, we study the effect of 
Ml and its phase ipu{i) on the decay widths of ti^2 and 61,2- In the numerical calculations, 
although the SUSY parameters Mi, M2, and Af are in general complex, we assume that 
ji, M2, At and Ajj are real, but Mi and its phase i^u{i) take values on the WMAP- allowed 
bands given in Ref. [5]. These bands also satisfy the EDM bounds [7]. And we evaluate 
the parameter M2 via the relation M2 = 3/5|Mi|(tan^vi^)~^. It is very important to insert 
the WMAP-allowed band in the plane Mi — ip into the numerical calculations instead of 
taking one fixed Mi value for all (/9-phases, because, for example, on the allowed band for 
= 200 GeV, Ml starts from 140 GeV for = and increasing monotonously it becomes 
165 GeV for (/? = tt. In Ref. [5] two WMAP-allowed band plots are given, one for ji — 200 
GeV and the other for ji — 350 GeV. For both plots the other parameters are fixed to be 
tan/3 = 10, mH+ = 1 TeV, At = 1.2 TeV, = 1.2 TeV, (^^=(/?^^=^^^=0. 
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2 Top and Bottom Squarks Masses, Mixing and De- 
cay Widths 

2.1 Masses and mixing in squark sector 

The superpartners of the SM fermions with left and right hehcity are the left and right 
sfermions. In the case of top squark (stop) and bottom squark (sbottom) the left and right 
states are in general mixed. Therefore, the sfermion mass terms of the Lagrangian are 
described in the basis {qL,qR) as 

2 :i ) C: ) 

with 

= M^ + {Ilj^-egSm^9w)cos{2p)ml + ml (2.2) 
Mrr = M^, +egSm'^9wcos{2(3)ml + ml (2.3) 
M'nL = (Ml^)* = m,(A,-/.*(tan/?)-2^3.) (2.4) 

where rUq, Cg, I^j^ and 9w are the mass, electric charge, weak isospin of the quark q=b,t and 
the weak mixing angle, respectively, tan/? = V2/V1 with Vi being the vacuum expectation 
values of the Higgs fields Hf, i = 1,2. The soft-breaking parameters Mq, Mq, = M^{Mf)) 
for q=t(b), Af, and At involved in Eqs. (2.2-2.4) can be evaluated for our numerical 
calculations using the following relations 

M? = ^ (ml + m| ± ^ (m| - m|)2 - 4mj\At - fi* cot /^p) 

-(^ - ^ sin" Ow) cos{2(3)ml - ml (2.5) 

Ml = ^ (ml + m| T ^ (m| - m|J2 - 4m|| A - cot /jp) 
2 

- - sin^ 9w cos{2(5)ml - mj (2.6) 

and similar ones for Mq and by interchanging t <-> 6 in Eqs. (2.5-2.6 ). 

The squark mass eigenstates gi and q2 can be obtained from the weak states and qu 
via the g-mixing matrix 

e*'^^ cos 9a sin 9a 
W=\ I (2.7) 

— sin 9q e ^'^^ cos 9q 



where 

if^ = arg[M2 J = arg[/l, - /,* {te^n (3)-'^ (2.8) 

and 




One can easily get the following squark mass eigenvalues by diagonalizing the mass matrix 
in Eq. (2.1): 

- I {mIl + ML T ^(M|^-M|J2 + 4|M|^|2^ , m^, < m,-, (2.10) 

We might add a comment about the possibility of a flavor mixing, for example, between 
the second and third squark families. In this case, the sfermion mass matrix in Eq. (2.1) 
becomes a 4x4 matrix in the basis (c^, cr, tL, tn). Then one obtains squark mass eigenstates 
(ci, C2, ii, from these weak states, and analyzes their decays by utilizing procedures 
similar to the ones indicated in the text. The problem with flavor violation effects is that 
their inclusion necessarily correlates B, D and K physics with direct sparticle searches at 
colliders. Moreover, it has been shown that, with sizeable supersymmetric flavor violation, 
even the Higgs phenomenology at the LHC correlates with that of the rare processes [8] . 
In this work we have neglected such effects; however, we emphasize that inclusion of such 
effects can give important information on mechanism that breaks supersymmetry via decay 
products of squarks. 

2.2 Fermionic decay widths of U and hi 

The quark-squark-chargino and quark-squark-neutralino Lagrangians have been first given 
in Ref. [1]. Here we use them in notations of Ref. [6]: 

= gt{l%PR + klPL)xpi + gWijPR + KjPL)xt% + h.c. (2.11) 

and 

= m{4kPR + biPL)xlqi + h.c. (2.12) 
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We also borrow the formulas for the partial decay widths of ( qi=ii and hi) into quark- 
chargino (or neutralino) from Ref. [6]: 



167rm| 



X 



m„ 



4i?e(4*4)m,.m,±]2.13) 



and 



g^X^/^ {ml, 7711,77110) 
167rm? 



X 



+ l4r)K -ml- m|o) - 4i?e(ag6f,)m,m,o] (2.14) 



with \{x, y, z) = + + z"^ — 2{xy + xz + yz). 

The explicit forms of fc^^ and a^^, h^^ can be found in Ref. [6]. We have to point 
out that although at the loop level the SUSY-QCD corrections could be important, our 
analysis here are merely at tree level, as can be seen from Eqs. (2.13) and (2.14). In this 
work we content with tree-level amplitudes as we aim at determining the phase-sensitivities 
of the decay rates, mainly. 



3 Top squark decays 

Here we present the dependences of the ii and ^2 partial decay widths on '{)u{i) foi' = 200 



GeV and for = 350 GeV. We also choose reasonable values for the masses (m^^, 



m 



t2l 



m^±. m-±, m^o 



(350 GeV, 800 GcV, 180 GeV, 336 GeV, 150 GeV) for 11 = 200 GcV 
and (mj^, m^^, m-±, m-±, m^o) = (350 GeV, 800 GeV, 340 GeV, 680 GeV, 290 GeV) for 
II = 350 GeV. 

For both sets of values by calculating the Mq and Mjy values corresponding to m^^ and 
mfj, we plot the decay widths for Mq > Mfj and Mq < Mfy, separately. Fig. 1(a) and 
Fig. 1(b) show the partial decay widths of the channels ii — > bxt , ii bxt h — > txi, 
h bxi, ^2 bxt and t2 tXi ^ function of (pu(i) for fi = 200 GeV assuming 
Mq > Mfj and Mq < Mjy, respectively. There we see some dependences on (pu{i) phase. 
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In order to see these dependences more pronouncedly, in Figs. 5(a)-(d) (ti 2 decays for 
H = 200 GeV ) and Figs. 6(a)-(d) (ti^2 decays for fi = 350 GcV ) we plot now separately 
only those decay channels whose dependences on fu{i) are not clearly seen in Fig.l and 
Fig. 2. r{ti txi) and r{ti bxf) decay widths increase as (pu{i) increases from to 
TT, but r(fi — > bxt) width decreases as (pu{i) increases. On the other hand, r(t2 — > tXi) 
decrease for both Mq > Mfj and Mq < Mfj cases as (pu{i) increases; r(i2 ^xt) decreases 
for Mq > Mjj but increases for Mq < Mjj and T{i2 — > 6x2 ) increases for Mq > Mfj but 
decreases for Mq < M^. 

The branching ratios for t2 are roughly B(t2 bxi) '■ -8(^2 ^X?) • B(t2 bx2) ~ 8 
: 2 : 1. This simply reflects both the large phase space and large Yukawa coupling for the 
decay ^2 ^ bxt- 

Figs. 2(a) and 2(b) show the same partial decay widths for fj, — 350 GeV. [See also 
Figs.6(a)-(d)]. They also reveal the significant dependences on CP- violation phase. 

For fjL — 350 GeV the WMAP-allowed band [5] takes place in larger Mi values (~ 
305 — 325 GeV) leading to larger chargino and neutralino masses. This naturally leads 
very small decay widths for ti — > bxt kinematically forbidden ti — * bxi and ti bxt . 
Because of the large t2 mass, ^2 decay processes are still kinematically allowed as seen in 
Figs. 2(a)-(b). The decay width of the process t2 bxt ^^c largest one among the t2 
channels and the branching ratios are B{i2 — > bxt) '■ -^(^2 bxt) '■ -^(^2 ^X?) ~ 4 : 2 
: 0.3. The decay ^2 — ^Xi shows strong phase dependence. 

4 Bottom squark decays 

We give sbottom decay widths as a function of i^u{i) in Figs. 3(a)-(b)(for ji — 200 GeV ) 
and in Figs. 4(a)-(b)(for ^ = 350 GeV ). [See also Figs. 7(a)-(d) and Figs. 8(a)-(c)]. Here 
we choose the masses (m^^, m^^, rn^±, m^±, m^o) = (550 GeV, 800 GeV, 180 GeV, 336 
GeV, 150 GeV) for fi = 200 GeV and (m^^, m^^, m-±, m-±, m^o) = (550 GeV, 800 GeV, 
340 GeV, 680 GeV, 290 GeV) for ^ = 350 GcV. 1(62 -> 6x?) is smaller than 1(62 ^ tx~) 
in spite of large phase space, because in 62 &Xi decay only Yi, coupling enters. The 
dependences of the phase (fiu{i) in sbottom decays are similar to those in stop decays. 
The branching ratios for 62 decays are B(b2 — > txi) ■ -B(&2 — ^ ^X2 ) ■ -^(^2 — ^'Xi) ~ 
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10 : 7 : 0.5. While the process 62 ^X? is suppressed more than one order its dependence 
on ipu{i) is prominent such that the value of decay width at <fu{i) = is 2 times larger 
than that at (pu{i) = tt. (^-dependences of the process 62 —>■ txi and 62 tX2 can be seen 
easily in Fig. 3(a). 

For MqKMf), bi — > tXi decay is about five times larger than 62 — tXi- The branching 
ratios for bi decays are B{bi txi) B{bi tx2) '■ B{bi bxi) ~ 5 : 1 : 0.2. 

5 Conclusions 

In this note, we have extended the study of third family squarks in MSSM with complex 
parameters in Ref. [6] taking into account the cosmological bounds imposed by WMAP 
data. For this purpose, we have calculated numerically the decay widths of the third family 
supcrpartncrs ti^2 and 61^2, in particular, their dependences on the CP phase ^pu{i)- We have 
found that some decay channels like ^2 — ^ bxt , ^2 i>X2 ^ ^2 tXi, ti — >■ bx^, bi — >■ tXi , 
bi bxi, 62 txi and 62 tX2 show considerable dependences on </?j7(i) phase. These 
decay modes will be observable at the LHC. Therefore they provide viable probes of CP 
violation beyond the simple CKM framework; moreover, they carry important information 
about the mechanism that brakes Supersymmetry. 

The HEP community eagerly waiting the discovery of supcrpartncrs at the LHC. The 
obvious task in this context is to recognize them, for example via their spins, experimen- 
tally. There has been several studies in this direction [9,10]. Recently, Wang and Yavin [11] 
have studied in detail the possibility of measuring the spin of new particles in a variety of 
cascade decays, especially in the decay channel q — > ?X^(x°) which we consider here, and 
they conclude that the prospects for spin determination are rather good. 
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Figure 1: (a)-(b) Partial decay widths T of the ti^2 decays for /i = 200 GeV , tan/3 = 10, 
At = 1.2 TeV, ip^=ipA=0, mi^ = 350 GeV and m^-^ = 800 GeV. 
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Figure 2: (a)-(b) Partial decay widths F of the t\^2 decays for /i = 350 GeV , tan/? = 10, 
At = 1.2 TeV, ip^=ipA=0 mi^ = 350 GeV and m^-^ = 800 GeV. 
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Figure 3: (a)-(b) Partial decay widths T of the 61^2 decays for fi = 200 GeV , tan/3 = 10, 
Ab = 1.2 TeV, ip^=ipA,=0 mi^ = 550 GeV and m-^^ = 800 GeV. 
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Figure 6: (a)-(d) dependences of certain t\^2 decays for \i = 350 GeV. 
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Figure 7: (a)-(d) (pu{i) dependences of certain 61^2 decays for /i = 200 GeV. 
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Figure 8: (a)-(d) fu{i) dependences of certain bi^2 decays for /i = 350 GeV. 
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